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Abstract 

We study a scenario that large non-Gaussianity arises from the baryon asymmetry of the Uni- 
verse. There are baryogenesis scenarios containing a hght scalar field, which may result in baryonic 
isocurvature perturbations with some amount of non-Gaussianity As an explicit example we con- 
sider the Affleck-Dine mechanism and show that a flat direction of the supersymmeteric standard 
model can generate large non-Gaussianity in the curvature perturbations, satisfying the observa- 
tional constraints on the baryonic isocurvature perturbations. The sign of a non-linearity parame- 
ter, /nl, is negative, if the Affleck-Dine mechanism accounts for the observed baryon asymmetry; 
otherwise it can be either positive or negative. 
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I. INTRODUCTION 

The WMAP results fl^ provided strong support for the inflation; the observed primordial 
fluctuations are consistent with nearly scale-invariant, adiabatic and Gaussian density per- 
turbations, as predicted by a simple class of inflation models. Those predictions are derived 
based on a simple but crude assumption that it is only the inflaton that acquires sizable 
quantum fluctuations during inflation. Its apparent success, however, does not necessarily 
mean that such a non-trivial condition is commonly met in the landscape of the inflation 
theory. 

It is perhaps natural to expect that there are many scalar flelds in Nature. If some of them 
are light during inflation, they acquire quantum fluctuations, which may leave slight devia- 
tion from the above properties in e.g. the cosmic microwave background (CMB) anisotropy, 
such as the isocurvature perturbations and/or sizable non-Gaussianity. Interestingly, it was 
recently reported that large non-Gaussianity was detected by the analysis on the WMAP 
3yr data [2J. The latest WMAP 5yr data seem to have the same tendency, although the 
vanishing non-Gaussianity is allowed within 95% C.L. [1]. Those hints on non-Gaussianity 
may be originated from such additional light scalars. 

The non-Gaussian features in the observed CMB can arise from either adiabatic or isocur- 
vature density perturbations. While the former was extensively studied in e.g. the curvaton 
0, 0]/ungaussiton [5] scenarios, much less attention was paid to the latter case since this 
possibility was noted in Refs. jo, Q]. 

Recently Sekiguchi, Suyama and the current authors systematically studied how the 
non-Gaussian isocurvature perturbations exhibit themselves in the CMB temperature fluc- 
tuations [8]. In particular, it turned out that the non-Gaussianity in the isocurvature per- 
turbations is enhanced at large scales, which may be conflrmed (or refuted) by the current 
and future observations. In Ref. [8] we also studied the QCD axion as an example, which 
generates isocurvature perturbations in the cold dark matter (CDM). In this paper we focus 
a scenario that baryonic isocurvature perturbations possess non-Gaussian properties. 

One of the promising candidates for providing a theory beyond the standard model (SM) 
is supersymmetry (SUSY), and the supersymmetric SM (SSM) contains many flat directions. 
Those flat directions can play important roles in cosmology; they generate the baryon asym- 
metry of the Universe via the Aflleck-Dine (AD) mechanism [9|, llO|, and may account for 
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dark matter by deforming into Q-balls jll, 12,[l^, [l^. If a flat direction remains light during 
inflation, it acquires quantum fluctuations, leading to the baryonic isocurvature fluctuations 
Il7l |. Moreover, as shown in [17], the phase direction of the flat direction generically 
remains flat in most inflation models in supergravity. Thus, we expect that the baryonic 
isocurvature perturbations as well as the associated non-Gaussianity are generically present 
in the AD mechanism. 

In this paper we study the non-Gaussian property of the baryonic isocurvature pertur- 
bations produced in the AD mechanism. We flnd that the resultant non-Gaussianity has 
distinctive features, compared to those produced in the curvaton/ungaussiton mechanism. 
In terms of a non-linerity parameter, /nl, to be deflned later, the AD mechanism predicts 
a negative value of /nl if the mechanism accounts for the observed baryon asymmetry of 
the Universe; otherwise, /nl becomes positive or negative, depending on the sign of baryon 
asymmetry created by the AD mechanism. We would like to emphasize that the AD mech- 
anism can generate a signiflcant amount of non-Gaussianity while accounting for the total 
baryon asymmetry of the Universe. 

This paper is organized as follows. In Sec. [Ill we briefly summarize the calculation of non- 
Gaussianity from isocurvature perturbations. Then we discuss non-Gausianity generated by 
the AD mechanism in Sec. [TTTl Sec. [IV] is devoted to discussion and conclusions. 



II. NON-GAUSSIANITY FROM ISOCURVATURE PERTURBATIONS 

We write the spacetime metric as 

ds^ = -M^dt^ + a^{t)e^^-i,j {dx' + f3'dt) {dx' + f3'dt) , (1) 

where M is the lapse function, (3i the shift vector, the spatial metric, a{t) the background 
scale factor, and the curvature perturbation. We denote by C the curvature perturbation 
ij) evaluated on the uniform-density slicing. The power spectrum PQ{k) and the bispectrum 
B(^{ki^ A;2, ks) of ( are deflned by the two-point and three-point correlation functions as 

(C^^C,^^) ^ {2nfS{h + k2)Pdh), (2) 
(C.-,C,-,C,-3) = i2nfSih + 4 + h) B^h, h), (3) 
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where is a Fourier component of i.e., C^. = J dPxe ^^^'^("(x), and ki = \ki\ for i = 1, 2, 3. 
The non-hnearity parameter /nl is defined by 

ks) = ^/NL[i"c(fci)n(^2) + 2 perms.]. (4) 

Let us now define the CDM and baryon isocurvature perturbations in the radiation- 
dominated universe as 

S,^ = 3{C,-C), (5) 
S,^ = 3{a-Q, (6) 

where (x is the curvature perturbation on a shcing where the energy density of the compo- 
nent X is spatiaUy uniform, and x = {c, 6, r} corresponds to CDM, baryon, and radiation, 
respectively. Since the baryonic isocurvature perturbation cannot be distinguished from the 
CDM isocurvature one, it is useful to define the effective CDM isocurvature perturbation as 

S = + RStj^ (7) 

where R Qi,/Qcdm — 0.2. 

We can similarly define the power spectrum Ps{k) and the bispectrum Bs{ki^ /c2, ks) of 
the eflFective CDM isocurvature perturbation S as 

{Sj^^Sj:^) ^ {2nfS{h + k2)Psiki), (8) 
{Sk,Sj:^Sj:^) = {27rf S{h + 4 + ^3) Bs{h, h), (9) 

and the non-linearity parameter fs is defined by 

Bs{kr, k2, ks) = fs[Ps{h)Ps{k2) + 2 perms.]. (10) 

Suppose that the isocurvature perturbation S is sourced by a scalar field, cj). Then S can 
be expanded in terms of the fiuctuation of (/) as 

where the fiuctuation Scj) is evaluated when the corresponding scale leaves the horizon during 
infiation. We define the power spectrum of the scalar field as 



= {27rfS{ki + k2)Ps^{ki). (12) 



In this paper we consider only the local type non-Gaussianity [1 



If the mass of (j) is much smaUer than i^inf, the power spectrum is approximately given by 

Psm ^ (13) 

where H[^f is the Hubble parameter during inflation, and we neglect the tilt of Hi^f for 
simplicity. For later use we also deflne the following: 



A?,.ilp«Wc.(|^y. (14) 

We can express Ps in terms of the Sep by substituting Eq. f fTTI ) into Eq. ([8]), 

Ps{k) ^ [si + Sl^A',^MkL)] Ps4k), (15) 

where we have introduced an infrared cutoff L, which is set to be of order of the present 
Hubble scale {3, Q, 2o|. In a similar way, we can express Bs in terms of Scj). The expression 
becomes simple when we take the so-called squeezed conflguration in which one of the three 
wavenumbers is much smaller than the other two (e.g. fci <C /c2, ^3), and it is given by 

Bs{h, k2, h) - [SlS^^ + Sl^/\]^ HhL)] [Ps^{h)Ps^{k2) + (2 perms.)] (16) 

where fc^ = min{/ci, /c2, /C3}. Thus, for the squeezed conflguration: ki <C ^2,^35 fs is given 

by B 

~ ^ fl7) 

- si + si,\Ai^\Hk,Ly ^ ^ 

In the following, we take the conflguration, ki <C ^2,^3, as the squeezed conflguration. 

Let us now relate fs to /nl- The curvature perturbation in the matter dominated era is 
given by 

C = C^^^ + (18) 

where ^^^^ denotes the primordial curvature perturbation created by the inflaton. We assume 
that the power spectrum of ( is predominantly produced by while the three-point 
correlation function originates from i.e., o:^ Bs/27. For the squeezed conflguration, 
fci <C /c2, fcs, we obtain 



Also) 5 ( Ps\^ r 
2 



162 V Pc 



C 



0. 



where we have used Eq. ( IT7I) . and we have defined = /c^P(^/(27r^) ^ 2.4 x 10 ^ [Ij. Here 
we have written the non-hnearity parameter as /^l^^ in order to emphasize that the non- 



Gaussianity comes from the isocurvature perturbation. This relation is insensitive to the 
wavenumbers, up to the tilt of the Pc^ and Ps- Note that the sign of /^^^ is determined by 
that of fs^ or equivalently, aS^^^^. 

Lastly we comment on the magnitude of /nl- The observed CMB temperature fiuctua- 
tions are consistent with the pure adiabatic perturbations, and there is a tiffht constraint on 
the isocurvature perturbations. According to the latest WMAP 5yr data the constraint 
reads Ps/ Pq^^-^^ at 95% C.L. for uncorrelated isocurvature perturbations. Therefore, in 
order to have large non-Gaussianity, |/nl| ^ 1, one can see from Eq. (fT9l) that \fs\ must be 
at least larger than 3 x 10^. It should be noticed however that the /nl in Eq. (fT9l) aflFects 
the CMB temperature fiuctuations in a completely different way from the conventional /nl 
defined for the adiabatic perturbation. That is to say, the currently available constraint on 
/nl, —9 < /nl < 111 at 95% C.L. [1], cannot be applied to the /nl^^ in our case. This is 
because the constraint is derived assuming that the non-Gaussianity arises from the adia- 
batic perturbations. What is more relevant to the CMB observations is /nl^, an eflFective 
non-linearity parameter defined by the three-point correlation function of the CMB tem- 
perature fiuctuations. For isocurvature perturbations with non-Gaussianity, sensitively 
depends on the scales of interest. Indeed, as pointed out in Ref. [s], f^^ is greatly enhanced 
as /nl^ ^ 100/|1Jl^^ large scales. No observational constraint on f^^ is known yet, and 
so, we estimate /nl*^ instead of /nl^, and take I/nl^^^I > 1 as the criterion for "large" non- 
Gaussianity. The reader should keep in mind that the f^^^ aflFects the CMB temperature 
fiuctuations differently from that defined for the adiabatic perturbations. Also note that 
although Eq. (fT8l) holds only for large scales in the matter dominated era, /nl°^ correctly 
characterizes small scale perturbations in the CMB anisotropy, once the transfer function of 
the isocurvature perturbation is taken into account, as explicitly shown in Ref. [st]. 



III. NON-GAUSSIANITY FROM AFFLECK-DINE MECHANISM 

There are several baryogenesis mechanisms containing a light scalar. Among them, we 
focus on the AD mechanism, which is described briefiy below. We will comment on other 
mechanisms in Sec. llVi 
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The SSM contains many flat directions consisting of squark, slepton and Higgs flelds. The 
flat directions are parameterized by composite gauge-invariant monomial operators such as 
udd or LHu^ and the dynamics of a flat direction can be expressed in terms of a complex 
scalar fleld 0, dubbed the AD fleld. The flat directions of the minimal SSM are classifled in 
Ref. [2j|. We assume that has a nonzero baryon number in the following. 



A flat direction has a vanishing scalar potential in the SUSY limit as long as only renor- 
malizable terms in the superpotential are considered, but it is lifted by a non-renormalizable 
operator in the superpotential: 

W^K = (20) 

where M denotes an effective cutoff scale for the interaction, and n is an integer: n = 4, 5, 6, 7 
and 9, which depends on flat directions. In addition to the non-reanormalizable operator, 
the flat direction is lifted by the SUSY breaking effects. 

In gravity-mediated SUSY breaking models, the scalar potential for the AD fleld (j) is 
given by 

Vs{<f.) = (ml - cH'M' + + h.c. j + (21) 

where c and are numerical constants of order unity, and is set to be real without loss 
of generality. Several comments are in order, m^^ is a soft SUSY breaking mass, and H is the 
Hubble parameter. We have included here the Hubble-induced mass, which arises from the 
quartic couplings between the AD fleld and the inflaton in the Kahler potential. We have 
assumed that the sign of the mass term is negative so that the AD fleld develops a large 
expectation value during inflation. Note that this mass term is present after inflation until 
the reheating is completed. The second term in Eq. f l2Tl ) is the baryon- number violating 
74-term, and the last one is due to the non-renormalizable operator f l20l ). We have dropped 
the so-called Hubble-induced 74-terms since they are absent in most inflation models in 
supergravity flT^. 



The AD fleld also feels flnite-temperature effects given by 



23| 



VTi<i>)= E Ckf,T'\<l>\' + aaiTfTHogO-^)^ (22) 

where is a constant of order unity, collectively denotes the gauge and Yukawa couplings 
for the corresponding AD field, and a is a constant of order unity. The sign of a depends on 
flat directions, and it is determined by the two-loop finite temperature effective potential. We 
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assume a to be positive in the following Those thermal effects are known to significantly 



aflFect the final baryon asymmetry 



26| 



Now let us take a closer look at the dynamics of the AD field. During infiation the AD 
field stays at the potential minimum, 

1^1 = (23) 

which is determined by the balance between the Hubble mass term and the non- 
renormalizable term. After infiation, the Hubble parameter decreases with time, and so 
does the minimum. The radial component of continues to track the minimum until it 
begins to oscillate. The oscillations start when the Hubble parameter becomes equal to i/os, 
given by 

^os = E ^kf.T^s + aaiTf^^. (24) 

Here the subscript "os" denotes that the variable should be evaluated at the beginning of 
the oscillations. The field value at which the AD field begins to oscillate is given by |(/)os| = 
(HosM^~^y/^^~'^\ On the other hand, the phase component of the AD field, 9 = arg[(/)], is 
almost massless during and after infiation, due to the absence of the Hubble-induced mass 
term. (Note that the mass along the phase direction arises only from the baryon-number 
violating yl-term.) Therefore 9 remains at the initial value set during infiation, until the AD 
field starts to oscillate. In the following we refer the initial value by 9. 

When the starts to oscillate, it is also kicked into the phase direction due to the baryon- 
number violating A-teim. Most of the baryon asymmetry is created at this moment. The 
angular momentum of the motion in the complex plane of (j) is related to the baryon number 
density: 

^ (25) 
The baryon-to-entropy ratio created by the AD mechanism is estimated as 

■sin(n^), (26) 



HIM 



If a is negative, the AD field may be trapped by the negative thermal logarithmicpotential. Some explicit 
examples of the flat directions having the negative corrections are given in Ref. [2J] . If the trap lasts long 
enough, the AD mechanism may not work [25|. 
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where A = |a^|(n — 2)/3(n — 3) is a constant of order unity. Note that 9 serves as a CP phase 
for the successful baryogenesis. Here we have assumed that the osciUations begin before the 
reheating completes. The AD mechanism can account for the observed baryon asymmetry 
of the Universe, ub/s 8.8 x 10~^^ [1], for appropriate choices of the cutoff scale M and 
the reheating temperature Tr. 

Now let us consider the fluctuations of the AD fleld [l^, 22]. Since the angular 
direction of the AD fleld does not receive sizable corrections during inflation [iT], it has an 
unsuppressed quantum fluctuation 56. We deflne the power spectrum of 50 as 

= {2Tifm + h)Pse{h). (27) 

The magnitude of the fluctuations is given by 

^'»-&'-^'^-{^^^- 

where |(/>inf| = {Hi^iM^~'^Y^^'^~'^\ This fluctuation of the AD fleld results in the baryonic 
isocurvature fluctuation through the AD baryogenesis mechanism, 



67 



71 ^ 

ncoi{nd)5d {59 f + 



(29) 



where r = /tib represents the fraction of the baryon number created by the AD mech- 
anism to the total baryon number. If the AD mechanism is responsible for the total baryon 
asymmetry, r equals to 1. In the presence of the other baryogenesis, r can take any values in 
principle, and it can be even negative. As we will see, the isocurvature perturbation and/or 
non-Gaussianity induced by the AD fleld can be signiflcant for both cases of r = 1 and 
|r| <C 1. 

The effective CDM isocurvature perturbation S = RS^j is similarly expanded as 

S = S050+^See{Sef + ---, (30) 



with 



Se = Rrncot{n9), (31) 
See = -Rrn^. (32) 



Eq. (flTj) can be rewritten as 

fs - 



^ [cot^ne) + n^\Al\ ln(A:2L))"' . (33) 



Using Eq. (fT^ . we obtain 

p(iso) 



„ > 2 

5 /A2 



fr ^Yk\^! {So + S'eel^hlHhL)) See, 



162 [-^J + n^lA^^I In(fciL)) . (34) 

Note that f^"°^ is negtive (positive) for a positive (negative) value of r. The constraint on the 
isocurvature perturbation, Ps/ 0.19 now reads 

[Rnrf {coi^ine) + n^lA^^I In(fcL)) A^^ < 4.6 x 10-^^ (35) 

In Fig. [1] we have plotted the contours of sgTi{r)f^^^ = —1, —10, —100 and —1000, together 
with the constraint on the isocurvature perturbations (l35l) . In order to avoid the constraint, 
cot(n0) and Aso are limited to the following ranges: 

^cot(ne) < O(O.l), (36) 

Ha,, < 0(10-^), (37) 



and the non-Gaussianity parameter f^^"^ is bounded as ^ 60. 

Lastly let us make a comment on a characteristic behavior of the baryonic isocurvature 
perturbations (l29l) . The second order term , which represents the non-Gaussian perturba- 
tion, dominates over the linear fluctuation in the limit cot(n^) ^ 0. In this limit the baryon 
number itself created by the AD mechanism approaches to the maximum (see fl26l)). In 
Fig. [2] we show schematically the baryon asymmetry generated by the AD mechanism as 
a function of 9. We have indicated the regions where /^l*^ is predicted to be positive or 
negative. This helps us understand why /^l*^ becomes negative when the positive baryon 
number is generated. 



A. Gravity-mediated SUSY breaking models 

In gravity-mediation models, the whole potential of the AD fleld is given by Eqs. (I2T1) 
and (l22l) . In the following we assume n = 6 For simplicity, we do not consider the case 
that oscillation due to the thermal mass term occurs. This assumption is justifled since 

In the case of n = 4, the allowed regions become smaller than those in the case of n = 6. 
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FIG. 1: Contour of sgn(r)/j 



(iso) 
NL 



-1,-10,-100 and -1000 (solid lines), together with the 



constraint on the isocurvature perturbation, Ps / Pc^^^.lQ (dotted line). The region above the 
constraint is excluded. We have set n = 6 and approximated ln(A:L) ^ 1. 



we are interested in the regime \(j)\ ^ T to avoid the gravitino overproduction [28|, l29| and 
hence particles that couple to the AD field cannot be thermalized. The AD field begins 
to oscillate due to thermal logarithmic term if the reheating temperature satisfies a 
condition (hereafter we take a = 1), 

1 /mWV^^ 



Tr > n = 



R 



a 



Ml ) 



1.1 X 10*^ GeV 



0.1 



a 



\3/4 / M \ 
1 TeV/ V 10^6 GeV y 



3/4 



(38) 



Otherwise it begins its oscillations by the soft mass term. Thus r = n^^^/n^ is estimated 



as 



0.27 



0.35 



1 TeV/ V 10^ GeV/ 



Tr WlTeV\ 



3/2 



0.1 



a 



m3/2 \ / lO^ GeV \ 
1 TeVy \ Tr ) 



M 



1016 GeV 



lO^^j ""^"^^ 

3 

sin ( 77,6*) 



for Tr<T^^ 



for Tr>T'^ 



(39) 
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FIG. 2: A schematic picture for the baryon asymmetry generated by the AD mechanism. We 
indicate the regions where f'^^^ is predicted to be positive or negative. Roughly speaking, the 
sign of is determined by the curvature of the curves, since the coefficient of the second order 
perturbations is determined by the second derivative of the baryon number. 



The CDM isocurvature perturbation S is calculated as 



S 



Vl TeV/ \ m^, ) 



3/4 



R 



W GeV 



1/2 



inf 



5 X 10-V3/4 ( — 
\ a 

,1/4 



for Tr<TI, 

1/2 



\V4 /lO^ GeV\ 



1 TeV 



1/4 



-R 



/ cot(ng) 
1014 GeVy V 10-2 

3/4 



if 



1014 GeV / 



cot (n6') 
10-2 



X sin^/*(n^) for Tr > T\ 



R 



(40) 

if the linear term in Eq. dominates. This is constrained as Ag < 2.2 x 10-^ from 
WMAP5 results, as stated before. Then from Eq. f|T^ we can calculate /jIjl"'' ^is 

,^ / cot(n^ ^^ ^ 
1014 GeVy V 10-2 



(iso) 
NL 



1 TeV/ V 104 GeV 



2^1TeV\V H-^^i Y /cot(n^) 



X sin2(n^) for Tr < Tf 



R 



-1 X lO-'r^ ( — 

a 



\ /lO^ GeV\ / H^^i V /cot(ni 



ITeV/l Tr Jv 1014 GeV j 1 10"^ 



X sin(n^) for Tr > T^ 



(41) 
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Thus depending on the sign of the CP phase (sin(n^)) (or equivalently the baryon asymme- 
try), the non-hnearity parameter can be either positive or negative. 

Here we mention the effect of Q-hall formation. It is known that if a scalar field has a 
conserved global U{1) charge and if the scalar potential becomes fiatter than a quadratic 
potential at larger field values, there exists a stable configuration of the scalar field, called Q- 
ball, whose stability is ensured by the U{1) symmetry [11]. In the context of AD mechanism, 
this conserved U{1) charge is the baryon number. In the AD mechanism, the scalar potential 
tends to be fiatter than the quadratic term due to the renormalization group effects. Thus 
Q-balls are generically formed in the AD mechanism {l2 . 

The charge of a Q-ball is given by 



1,Q,[l4l. 



Q - 7 — 



X 



for e > Cc 
for e < Cr 



(42) 



where 7 ^ 6 x 10~^ and Cc ^ 0.01 from the lattice calculation [311. The parameter e is called 
the ellipticity parameter, defined by the ratio of the minor and major axes of the orbit of 
the AD field (/). It is roughly estimated as 



ms/2 . . . 
6 ~ sm[nu). 



(43) 



If (5 ^ 10^°, the decay temperature of the Q-ball becomes smaller than the freeze-out 
temperature of the lightest supersymmetric particle (LSP) and hence LSPs emitted by the 



Q-ha 
LSP 



1 d ecay m ay b e overproduced, depending on the self-annihilation cross section of the 



jcav 

14, y, 



3^ 



Q 



3j]. The charge of the Q-ball, is estimated as 
/ \ / 1 Te V \ ^ / ]\/[ \ ^/^ 



2 X 10^1 



O.n^ /lO^ Gevy 



for Tr<T] 



-R 



M 



(44) 



for Tu>T: 



-R 



Tr J GeV^ 

If some fraction of LSP dark matter comes from the Q-ball decay, it produces the additional 
CDM isocurvature perturbation with some amount of non-Gaussianity. Furthermore, the 
decay rate of the Q-ball depends on its charge Q and hence the decay rate also becomes 
the source for the isocurvature fiuctuations [3&\ as well as non-Gaussianity. However, in 



Precisely speaking, the U{1) symmetry is explicitly broken by the A-term in the AD mechanism. Never- 
theless the effect of the breaking is very small after the AD field starts oscillating [30] . 
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Tr [GeV] 




Hint [GeV] 

FIG. 3: Constraints on {H^^i^Tji) for 7713/2 = 1 TeV and coi{n9) = 10~^. We have fixed r = 1. 
The lines show : f^^^ — —1 (red sohd), —10 (purple solid), isocurvature bound (green dashed), 
gravitino overproduction bound (blue dotted), LSP overproduction from Q-ball decay (orange 
long-dashed). Arrows indicate the allowed region. 

the interesting parameter regions, we have checked that the Q-balls evaporate in the high- 
temperature plasma and hence those effects of the Q-ball decay are negligible. 

In Fig. [3] the contours of f^^^ = — 1 and —10 are shown by the red solid and purple solid 
lines on the (i/mf , Tr) plane for 7773/2 = 1 TeV and cot (77^) = 10~^. Here we have fixed r = 1, 
i.e., the AD mechanism creates the total baryon asymmetry of the Universe. We also show 
the constraints from the isocurvature perturbation (green dashed), baryon overproduction 
(brown dot-dashed), gravitino overproduction (blue dotted), and LSP overproduction from 
Q-balls (orange long-dashed). To be conservative, we have estimated LSP abundance by 
neglecting the annihilation of the LSP, and we set the LSP mass to be 100 GeV. One 
can see that a significant amount of non-Gaussianity can be generated without conflicting 
the isocurvature constraint for the Hubble scale during inflation H-^^i > 10^^ GeV. For 
cot (77^) ^0.1, the isocurvature constraint becomes more stringent and large /^l^^ cannot be 
generated without conflicting the isocurvature bound (see also Fig. [1]). 

In the above arguments we have assumed that the AD mechanism provides the total 
baryon number of the Universe. However, the AD mechanism may create only small fraction 
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of the baryon asymmetry, most of which is dominantly generated by another baryogenesis 
mechanism. In this case r can be much smaUer than unity. In Fig. H] we similarly show the 
constraints on (Hi^f^Tji) for r = ±10~^, 7713/2 = 1 TeV and cot{n9) = 10~^. If r is positive 
(negative), /^^^ becomes negative (positive). Thus large positive value of f^^\'^ 1) can be 
obtained if the AD mechanism creates small amount of baryon asymmetry with the negative 
sign. In this case the larger inflationary scale (i/inf ^ 10-^^ GeV) than the case of r = 1 is 
required (see Eq. fl34l)). 

It may also be useful to rewrite Eqs. f l40l ) and f l4TI ) by substituting f l39l ). 
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for rR<r^ 
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X sin(n^) 



O.iy ^ m3/2 \ / lO^ GeV \ / M \ i jij^f 

ITevH yVlO^'^GeVy V 10^^ GeV 



9/4 



3/4 



cot(n^) 
10-2 



for rR>r^ 



(45) 



and 
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(46) 

We also show the constraints on (M, Tr) plane in Fig. [5] without flxing r. Here we take 
m3/2 = 1 TeV, cot{ne) = lO'^ and = 10^^ GeV. The red sohd line shows fj-'^^ = -1. 
One can see that, for |r| < 1, the cutoff scale M should be smaller than around 10^^ GeV 
for generating large non-Gaussianity. 



B. Gauge- mediated SUSY breaking models 

In gauge-mediated SUSY breaking models (GMSB) {36^, messenger flelds mediate the 
SUSY breaking effect to the SSM sector. Here we assume the direct mediation scenario, 
that is, the SUSY breaking fleld S couples to messenger flelds (^^, ^) in the superpotential 
as tU = S'^^^^. In this type of models, the gauge mediation is suppressed for the scale |(/)| > 
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FIG. 4: Same as Fig. [HJ but for r = ±10 ^. The sign of /11 depends on the sign of r 



{S) where {S) is the vacuum expectation value of S. Instead there appears a logarithmic 



correction to the potential as 



37] 



={ml/,-cH'M' + V, (log IP)' 



+ am"^3/2 



nM 



n— 3 
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|0|2(n-l) 
M2("-3) 



(47) 



(1^1 > {S)), 



where Vq is given by 



(48) 



where ttj is the gauge couphngs relevant for the AD field with i denoting the gauge groups 
and Fs = V^m^/2Mp for vanishing cosmological constant. Thus the whole potential of the 
AD field is the sum of Eq. (1221) and (jlTI) . The AD field begins to oscillate due to thermal 
logarithmic term for 

(49) 
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FIG. 5: Constraints on [M.Tr) for mg/a = 1 TeV, i?inf = 10" GeV and cot(n6') = 10"^ The 
lines show : — —\ (red solid), isocurvature bound (green dashed), gravitino overproduction 

bound (blue dotted), baryon overproduction (brown dot-dashed) and LSP overproduction from 
Q-ball decay (orange long-dashed). Arrows indicate the allowed region. 



For Tr < T^, the AD field begins to oscillate due to the logarithmic potential coming from 
the gauge mediation effects.*^ The resultant baryon asymmetry can be estimated as 

m3/2 ) \W 



1.3 X 10"^ 



for Tr<T^ 



a 



.1 MeV. 

Thus the isocurvature perturbation is calculated as 
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3/4 



(50) 



(51) 



We found that in the interesting parameter regions in the following, the oscillation due to the 7713^2 1 
term does not occur. 
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if the linear term in Eq. (I29l) dominates. The non-hnearity parameter can be calculated as 
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(52) 



If the parameters take appropriate values, /^l^^ can be large enough with either positive or 
negative sign. 

In GMSB models, there also exist a Q-ball solution, but its properties differ from those 
of the gravity-mediation type. The charge of the Q-ball is estimated as 



Q ^ P 



Mf 
0.06 



for e > Q 
for e < er 



(53) 



where /3 ^ 6 x 10""^ and Cc ^ 0.06 [3j| and Mp in the denominator should be replaced 
with Tos if the oscillation begins due to thermal effects. If the Q-ball energy per unit charge 
{MqIQ ~ MfQ~^^^) is smaller than the nucleon mass ^ 1 GeV, Q-balls become stable 
against the decay into nucleons and contribute to some fraction of the dark matter |38l |. 
Also in that case only the evaporated charges in the high-temperature plasma provides the 
baryon number remaining in the universe 0, . In the regime of e < 6c the Q-ball charge 
can be calculated as 

m3/2 ) VlO'*' GeVy 

for Tr>T\ 
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5 X 10^^ 
2 X lOii 



0.1 

a 



-R 



^ ' 10^ GeV y 
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(54) 



One can check that 1 GeV < Mq/Q is always met in the parameter regions of our concern. 
Thus Q-balls are unstable in the most of the interesting parameter region and Q-ball forma- 
tion does not affect following results. (Note that no LSPs are produced by the Q-ball decay 
since it is kinematically forbidden.) 

Assuming that the AD mechanism creates the total baryon number of the Universe (r = 
1), we obtain constraints on (i/mf, Tr) plane as shown in Fig. [61 We have set 7713/2 = 10 GeV 
and cot(n^) = 10~^. Meanings of lines are same as Fig. [31 We can see that a significant 
amount of non-Gaussianity can be generated, similar to the gravity-mediation case in the 
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FIG. 6: Same as Fig. [2l but for the gauge mediation. We have set — 10 GeV, r — 1 and 

cot(n6') = 10-2. 



previous subsection. However, as we decrease the gravitino mass, the allowed region wiU 
become smaller. 

We are also interested in the case that the AD mechanism generates a small fraction 
of the baryon asymmetry, with both positive and negative signs. Substituting (l50l) into 
Eqs. dsn and (IS2D yields 
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FIG. 7: Same as Fig.[5l but for the gauge mediation. We have set 7713/2 = 10 GeV, i^inf = 10^^ GeV 
and cot (n(9) = 10"^. 

The resultant constraints are shown in Fig. [71 on [M^Tr] plane. We have set 7773/2 = 
10 GeV, i^inf = 10^^ GeV and cot(770) = 10~^. We can see that either positive or neg- 
ative /nl^^ can be obtained through the AD mechanism while satisfying the isocurvature 
constraint. 



IV. DISCUSSION AND CONCLUSIONS 

We have focused on the AD mechanism in this paper, but there are other baryogenesis 
scenarios that may induce large non-Gaussianity in the baryon asymmetry. For instance, 
it is known that the spontaneous baryogenesis j4lt] can lead to the baryonic isocurvature 
perturbations, since in the original model a non-vanishing chemical potential is due to a 
slow-rolling scalar field. In order to estimate non-Gaussianity, however, one has to specify 
how the chemical potential for the baryon number arises. In the case of the spontaneous 



baryogenesis using a fiat direction of the SSM [421], the resultant non-Gaussianity has almost 



the same features as that in the case of the AD mechanism. Another example is non-thermal 
leptogenesis using the right-handed sneutrino condensate, N Q]. Suppose that N is light 
and fiuctuating around the origin during infiation. Such a set-up may occur without fine- 
tuning, because the origin is the symmetry-enhanced point, and this is exactly what is 
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considered in the ungaussiton scenario [5^. The N can generate non-Gaussianity in both 
the adiabatic and baryonic isocurvature perturbations. If the mass of the right-handed 
sneutrino is heavy, the baryonic isocurvature perturbations wiU become more important 
than the adiabatic one. In this scenario, the non-Gaussianity becomes large only if the 
non-thermal leptogenesis accounts for a tiny fraction of the baryon asymmetry. 

It is also possible to consider non-Gaussianity in other type of isocurvature perturbations. 
For instance, a large lepton asymmetry may have isocurvature perturbations with some 



amount of non-Gaussianity. If they are generated by the AD field with a lepton number |43l |. 
non-Gaussianity arises from the fiuctuations of the phase component, as in the case of the 
AD mechanism. However, such isocurvature perturbations will aflFect the CMB temperature 
fiuctuations in a different way, and so does the associated non-Gaussianity. We leave this 
issue for future work. 

So far we have considered up to the three-point correlation functions. It is straightforward 
to extend our analysis to the correlation functions of higher order. In particular, when the 
linear perturbation is negligible, we have a consistency relation between /nl and tnl as 
in the case of the ungaussiton [5] , the latter of which is defined as a non-linearity parameter 
of the four-point correlation function. 

In summary, in this paper we have studied a scenario that non-Gaussian baryonic isocur- 
vature fiuctuation is produced from the AD mechanism in SUSY. We have found that the 
AD mechanism can create large non-Gaussianity without confiicting with the current isocur- 
vature constraint. We have seen that for I/nl^^^I ^ 1 the Hubble parameter during infiation 
must be larger than about 10-^^ GeV, and that is bounded as < 60 to satisfy 

the isocurvature constraint. Interestingly, as opposed to many known mechanisms for gen- 
erating sizable non-Gaussianity the AD field can generate large non-Gaussianity even if 
it is the main component for generating the baryon asymmetry. The non-linerity parameter 
/nT^ is negative if the AD mechanism is responsible for the total baryon asymmetry of the 
Universe, while it can be either positive or negative otherwise. 

In our scenario the non-Gaussianity is necessarily accompanied with some amount of 



To generate large non-Gaussianity while satisfying the constraint on the amplitude of the isocurvature 
perturbation, the linear term in Eq. (pTj) must be suppressed (see also Fig. [TJ) 

It is pointed out that negative /nl can be obtained in the curvaton scenario if the curvaton potential 
deviates from quadratic one 
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the isocurvature perturbations, therefore, if there are indeed large non-Gaussianity arising 
from the isocurvature perturbations, the future (or even on-going) observations wiU detect 
the isocurvature perturbations. It would be very interesting if the non-Gaussianity tells us 
about the origin of the baryon asymmetry, which is difficult to be probed, otherwise. 
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